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Abstract. Self-gravitating accretion discs have only been studied in a few nearby objects using maser spots at the parsec-scale. 
1 We find a new spectral window for observing the self-gravitating accretion disc in super-Eddington accreting Active Galactic 

Nuclei (AGNs). This window is determined by calculating the outermost radius (r sg ) of a non self-gravitating disc and the 
corresponding emission wavelength (A sg ) as a function of various disc parameters. We find that A sl! reaches ~ 4000A for a = 0.1, 
when M > 70(Mbh/10 7 M g F' L^m/c 2 (where a, M, Mbh and Lgdci are, respectively, the viscosity parameter, gas accretion rate 
onto the central black hole (BH), the BH mass and the Eddington luminosity). Moreover, A sg is as small as ~ 1500A for 
a = 0.001, which is the smallest a case in this study. Therefore, the window for observing the self-gravitating part of an AGN 



CN 



accretion disc is from ~ 2/im to A sg . Incidentally, r sg can be less than the photon trapping radius for M > 10 3 3 LEdd/c 2 . Namely, 



'sg < 

a self-gravitating, optically-thick, advection-dominated accretion disc is expected to appear in the extremely high accretion rate 
regime. 

' Next, we demonstrate that the Mid-Infrared to X-ray spectrum of a bright, well-studied Narrow-Line Seyfert 1 galaxy, 

Ton S 180, is indeed well fitted by the spectrum arising from the following three components: an inner slim disc (with a 
corona), an outer, self-gravitating non-Keplerian disc and a dusty torus. The total mass, BH mass plus the entire disc mass, is 
found to be about (1.4 - 8.0)Mbh- If the surface density varies with radius r in proportion to r~ 06 , the total mass is consistent 
with the central mass estimated by FT/? and [O III] widths. 

o , 
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1. Introduction & Boiler 1998; Mineshige et al. 2000). The gas accretion 

«_J ■ rate M for those objects is expected to be comparable to or 

^ The Optical/UV/X-ray emission from Active Galactic Nuclei , arger to 10LEdd/c 2 5 where ^ is me Eddington luminos . 

(AGNs) is thought to arise from a hot accretion disc around ky For such super . Eddington accret ion rates, the outer edge of 

a super-massive Black Hole (BH; e.g., Mushotzky, Done & ^ nQn self . gravitating part of the accret ion disc radiates opti- 

Pounds 1993; Koratkar & Blaes 1999), while the Mid- and cal continuum emission (~ 4000A) for a certain parameter set 

Near Infrared (IR) radiation comes from a cold, dusty torus (Kawaguchi 2003) . they ^ p 0teri tially good candidates 

(e.g., Telesco et al. 1984; Thatte et al. 1997; Pier & Krolik tQ smdy ±e self . gravitating part of the disc obs ervationally at 

1993). Disc self-gravity can generally be ignored in the inner- longer wavelengths (Collin et aL 2002; Kawaguchi 2003). 
most region where the high energy spectrum is formed. At large 
distances (typically > 10 3 — 10 4 Schwarzschild radii, depend- 
ing on accretion parameters), the disc mass is expected to play 

a significant role. Disc self-gravity has been tested only at the In this paper, we qualitatively describe the wavelength cor- 
parsec-scale via maser spots detected in a few nearby objects responding to the emission from the outer edge of a non self- 
(e.g. Nakai, Inoue & Miyoshi 1993; Miyoshi et al. 1995; Hure gravitating disc for various accretion parameter sets. This is 
2002; Lodato & Bertin 2003). done in §2. We present in §3 a Mid-IR to X-ray spectral mod- 
Narrow Line Seyfert 1 galaxies (NLSls) and their high lu- eling of a bright, well-studied NLS1, Ton S 180 (PHL912), 
minosity analogue, Narrow-Line QSOs, are supposed to have for which numerous, multi-waveband observations are avail- 
high accretion rates among the AGN population (e.g, Brandt able (e.g., Wisotzki et al. 1995; Turner et al. 2002; Vaughan 
et al. 2002). In §4, several discussions are presented. The final 



Send offprint requests to: T. Kawaguchi section is devoted to a summary. 



2 



T. Kawaguchi et al.: Optical emission from super-Eddington accreting AGN 
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non self-gravitating, 
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self-gravitating, thin disc ou + 




Fig. 1. Schematic configuration (not to scale) of the three rel- 
evant components surrounding the black hole: a slim disc, a 
self-gravitating disc and a torus. For Ton S 180 (§3), the char- 
acteristic radii are found to be the followings: r sg 3 x 10 3 rs c h, 
font - 3 x 10 4 rsch and r x 3 x 10 5 rs c h- The broad-line region is 
located at r 1.4 x 10 5 r Sc h (§4). 

2. Critical Radius and Critical Wavelength 

In this section, we calculate the radius of the outer edge of the 
non self-gravitating slim disc (r sg ), i.e. inner edge of the self- 
gravitating disc, as a function of the BH mass Mbh, M and 
viscosity parameter a. We will then show the corresponding 
wavelength (A sg ) that is relevant to emission from r sg . We define 
M 7 for M B h in the unit of 1O 7 M . 

Figure ^ illustrates the configuration considered in this 
study. 

2.1. Critical Radius 

For super-Eddington accretion rates (i.e. M > 10LEdd/c 2 ), 
an optically-thick, advection-dominated regime ap- 
pears (Abramowicz et al. 1988). Inside a certain radius 
r trap ~ 0.5M(LEdd/c 2 ) _1 r S ch (Begelman & Meier 1982; 
Kawaguchi 2003), the timescale for photon diffusion is larger 
than the accretion timescale, and thus most of the photons 
emitted inside the flow are trapped within it. Thus, advection 
cooling becomes dominant over radiative cooling at r < r tlap . 

We define r sg as the radius where the mass density in the 
disc mid-plane (p m id- P iane) equals p sg . Here, p sg is a critical den- 
sity above which self-gravity of the disc must be considered: 



p sg = Q.l/(4nG) = 2.2M 7 - 2 (r/3r Sch )- 3 g/cm 3 , 



(1) 



where Qk is the Keplerian rotation frequency (e.g., Goldreich 
& Lynden-Bell 1965; Hure 1998). Inside r sg , Pmid-piane i s 
smaller than p sg , thus the disc self-gravity is negligible there. 

Many authors have tried to derive r sg (e.g., Shore & White 
1982). Actually, r sg (MBH, M, a) derived and used in this study 
is qualitatively consistent with that in Hure (1998). 

In order to derive r sg , we use analytical formulae for the 
mass density (p m id- P iane) given for a Newtonian gravitational 
potential (e.g., Kato, Fukue & Mineshige 1998). We confirm 
that these standard-disc formulae for various physical quanti- 
ties are consistent with numerical computations (Kawaguchi 
2003) within a factor of ~ 1.5 at any radius larger than r^. 
In other words, the flow behaves like a standard disc (Shakura 
& Sunyaev 1973) outside r^, even for a super-Eddington ac- 
cretion rate. 
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Fig. 2. Classification of the region where r sg appears (§2.1). 
Region (c) harbors r sg when M is less than M c b, while re- 
gions (b) and (a) are, respectively, in charge of r sg when M 
is smaller and larger than Mba- Around the radius r trap , the flow 
undergoes a change from radiation cooling-dominated regime 
to advection-dominated regime (optically-thick ADAF). 



In a standard disc, the following three regions must be dis- 
tinguished according to the sources of the pressure and opacity 
(e.g., Shakura & Sunyaev 1973; Kato et al. 1998). In which re- 
gion r sg appears is determined by M (see Figure|2)- 

- an outer region (c), where the gas pressure and absorp- 
tion opacity dominate over radiation pressure and electron- 
scattering opacity, respectively. With a relatively low ac- 
cretion rate, i.e. M < M cb = 2M^ 1 (a/0.1) 7/13 iEdd/c 2 , r sg 
appears at the outer region (c): 



'sg,c 



\0.lj 7 U Ed d/c 2 / 



-22/45 



rsch- (2) 



- a middle region (b) where gas pressure and electron- 
scattering opacity are dominant. For a moderately high ac- 
cretion rate, i.e. M < M ba = 70M- 1 (a/0.1) a4 L E dd/c 2 , r sg 
is given by: 



,14/27 / x -8/27 

" \ ■ T -26/27 / M 



,,,. = 5700 -1 ^"{J^p 



rsch, 



(3) 



- an inner region (a) which is described by radiation pressure 
and electron-scattering opacity. For a higher accretion rate 
(M > Mba), 7"sg is described as 



, sg , a = 250 - 



2/9 



2/9 / M 



,L Edd /c z 



4/9 



rsch- 



(4) 
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Given M B h and M, the critical radius r sg can be written 
as Max(r sg c , r sgj b, ?"sg,aX f° r the whole range of accretion rates. 
Instead of a piece-wise function, we employ a simple formula: 
r sg - ( r sg,a 3 + r sg,b 3 + '"sg.c 3 ) 1 ' 3 - Figure |3 (upper panels) shows 
r sg for various M BH , M and a. Until the accretion rate reaches 
Mba, r sg decreases with an increasing M, while r sg changes the 
sign of its M-dependency at M ~ Mba- 

We find that with M > 10 33 (a/0.001) 04 
(M B h/1O 9 M ) -0 - 4 Lndd/c 2 , r sg can be less than r tmp . The 
above expressions are no longer valid in the regime in which 
'sg < ''trap- In that case, a self-gravitating, optic ally- thick, 
advection-dominated accretion disc appears. Therefore, we 
restrict ourselves to accretion rates lower than 10 3 5 LEdd/c 2 , in 
order to ensure that r S g > r^p. 



2.2. Critical Wavelength 

Next, we evaluate the effective temperature at r sg , and subse- 
quently a critical wavelength A sg corresponding to the emission 
from r sg . At r > r trap , the effective temperature 7^ at radius r 
is given by 



T e «(r) = 6.2 x 10 5 M: 



■1/4 



M 



^Edd/C 2 



1/4 



- 

\>"Sch 



-3/4 



K. 



(5) 



The simplest procedure to estimate A sg is to use the Wien 
displacement law assuming that the emission is close to a 
blackbody radiation: namely A sg =s hc/4kT e g (see Rybicki & 
Lightman 1979). It follows that the critical wavelength is 



he 



4£r eff (r sg ) 



0.36 



K 



cm. 



(6) 



The results are shown by the dotted lines in Fig.[3](lower pan- 
els). 

The effects of electron scattering (opacity and 
Comptonization) must be considered in studying the emergent 
spectra from a super-Eddington accreting disc (Shimura & 
Manmoto 2003; Wang & Netzer 2003; Kawaguchi 2003). 
Comptonization is important at radius less than ~ 300rs c h 
for M < 1000LEdd/c 2 , and thus it has no effect on the 
emission beyond r sg . On the other hand, the ratio of electron- 
scattering opacity to absorption opacity is larger than unity at 
r < 2500[M/(L E dd/c 2 )] 2/3 (see Kato et al. 1998), and thus the 
effect of electron-scattering opacity is still important at r sg . 

Due to the electron-scattering opacity, the emergent spec- 
trum (I v ) is distorted towards higher energy (i.e. shorter wave- 
length), and is close to the so called modified blackbody 
(Rybicki & Lightman 1979; Czerny & Elvis 1987; Wandel & 
Petrosian 1988): 

2 By(T e ff) 



1 + V^abs.v + K es )/Ka- 



(7) 



where B v , Afabs.v and /e es (= 0.4cm 2 /g) are, respectively, the 
Planck function, absorption opacity and electron-scattering 
opacity. The spectral shift due to this effect is roughly estimated 
as follows (Madej 1974): 

^abs(f" sg) 



A log 



10 ^sg 



- I logi 



Kabsfog) + K es (r sg ) 



(8) 



where A- a b s is the Rosseland mean opacity for absorption 
(bound-bound, bound-free and free-free processes). In assess- 
ing K a b s , we take 30 times the free-free absorption opacity, as is 
done in Czerny & Elvis (1987) and Kawaguchi (2003): 



K abs = 30/Cfr 

with the free-free absorption opacity, Kg, as follows 

22 / Prnid-plane \ ( ^mid-pla 



Kg = 6.4 x W 



-3.5 

ilane \ 2 



K 



cm /g. 



(9) 



(10) 



Here, T^d-piane is the temperature at the mid-plane of the disc 
(Kato et al. 1998). 

Finally, A sg modified from Eq. (|6j by Eq. are plotted in 
the lower panels of figure|5]as solid lines. As far as M < Mb a , 
A sg is roughly proportional to M~ - 5 . ForM > Mba anda = 0.1, 
A S g reaches and stays at ~ 4000A, and interestingly it can be as 
small as ~ 1500A for the lowest a case here (i.e. a = 0.001). 
Because of dust sublimation at T > 1500 K, emission of 
dust falls within the Mid-IR-Near-IR spectral range and can 
not contribute to the optical component (see §3.3). Therefore, 
emission from ~ 2/im to A sg arises from the self-gravitating 
part of the super-Eddington accreting disc. In other words, this 
is a discovery of a spectral window for observing the self- 
gravitating disc, which has solely been studied by maser spots. 

We note that A sg for Mbh S 10 9 M o is also small even with 
sub-Eddington accretion rates. In the sub-Eddington regime, 
however, heating by central radiation onto an outer region at 
r ~ r sg may not be negligible. An uncertainty of physical 
quantities (such as the gas density) in the outer region of 
the disc (e.g., H5shi & Inoue 1988) will lead to uncertainties 
of r sg and A sg . On the contrary, a super-Eddington accreting 
disc is not flared, i.e. it has a maximum of the aspect ratio at 
r ~ (5 - 100) rsch (Kawaguchi 2003), when such irradiation is 
considered to be less important. For this reason, r sg and A sg ob- 
tained for high M are more reliable than for the sub-Eddington 
cases. Therefore, we concentrate on the super-Eddington cases. 

3. Spectral Calculations 

We now demonstrate that the broad-band spectrum of a bright, 
well-studied NLS1, Ton S 180, is indeed well fitted by the 
spectrum arising from three components: a slim disc (with a 
coronal, hard X-ray spectrum), a self-gravitating non-Keplerian 
disc and a dusty torus. The configuration we consider is 
schematically shown in Fig.[0 

The Near-IR to X-ray data are taken from Fig. 7 and Table 6 
of Turner et al. (2002). The IRAS data for Far-IR and Mid- 
IR flux are derived from the NED database. The luminosity 
is calculated from observed flux assuming isotropic radiation, 
zero cosmological constant, deceleration parameter qo = 0.5, 
and Hubble constant Hq = 75 km/s/Mpc. Data are plotted in 
Fig-EJas open squares. Ton S 180 has low Galactic and intrinsic 
extinction (Turner et al. 2002). Since Ton S 180 is a nearby 
object (z=0.06), we have not taken any spectral shifts or K- 
corrections into account. 

The BH mass estimated from the HB width and optical lu- 
minosity (Wandel, Peterson & Malkan 1999; Kaspi et al. 2000) 
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1 10 100 1000 1 10 100 1000 1 10 100 1000 

M [L Edd /c*] M [L^/c*] M [L Edd /c*] 

Fig. 3. Critical radius and wavelength [r sg (upper) and A sg (lower)] for various M and a. Here, (dashed lines in the upper 
panels) is the radius inside which photon trapping plays a role. The dotted lines in the lower panel are computed based on the 
assumption that the radiation from r sg is blackbody, while solid lines take the effect of electron scattering (modified blackbody) 
into account. If A sg < 2/j.m, it means that we are able to observe the self-gravitating accretion disc as Near-IR/optical continuum 
radiation. 



is 10 7 l M Q . This value is comparable to 10 73 M Q inferred from 
the [O III] width using Eq. (2) in Wang & Lu (2001). Its B- 
band luminosity, vL v (B), is about 10 446 erg/s, implying that 
M ^ 500LEdd/c 2 (see Kawaguchi 2003). Thus, this object is 
expected to be one of the highest M /(Lndd/c 2 ) objects among 
NLSls and Narrow-Line QSOs. 

3.1. Spectrum of the inner, slim disc 

The structure and emergent spectrum of the inner disc around a 
non-rotating BH is computed numerically following the same 
procedure as in Kawaguchi (2003). The theoretical background 
of the numerical method is described in Matsumoto et al. 
(1984) and in Honma et al. (1991a,b). The model employs 
the standard prescription for turbulent viscosity v t = ac s H 
(Shakura & Sunyaev 1973) where c s is the sound speed and 
H the semi-thickness. The input parameters are Mbh, M and 
viscosity parameter a. Given a set of these parameters, we in- 
tegrate the differential equations for the flow from 2 x 10 4 rs c h 
down to the inner free boundary at 1.01rs c h- 

The best parameters for this object have been selected by 
the least-square procedure, comparing the predicted spectrum 
(with a hard X-ray power-law due to a corona) to the power- 
law fits for the optical/UV/X-ray spectrum (thick striped lines 
in Fig.|4j see Table 7 in Turner et al. 2002). 

First, the following parameter region is surveyed; Mbh = 
lO 6 - 5 - 7 - 5 M , M = 100 - 1000L Edd /c 2 , and a = 0.1 - 0.001 



Table 1. Best parameter sets for the non self-gravitating disc. 



set 


M BH [M G ] 


rsch[cm] 


M[L Edd /c 2 ] 


r S g[>sch] 


£[g/cm 2 f 


#[cm] b 


1 


10 7 


3.0 10 12 


10 2.75 


2000 


2.4 10 6 


1.2 10 15 


2 


10 6.8 


1.9 10 12 


1000 


3000 


2.6 10 6 


1.3 10 15 



a Half surface density £ at r = r\ 
b Semi thickness H at r = r s „. 



with Alog(M B H) = Alog(M) = 0.5, and Alog(o-) = 1. Among 
them, one parameter set (M B h = 10 7 M o , M = 10 2 5 L E dd/c 2 and 
a - 0.001) exhibits the best fit. Next, we compute the spectra in 
more detail around this parameter set, with Alog(MBH) = 0.2, 
Alog(M) = 0.25, and Alog(o') = 0.3. Finally, we selected two 
best sets of parameters as listed in Tabled with a = 0.002 in 
both cases. The mass of this inner disc Mdi sc (r < r sg ) is 0.009 x 
M B h for set 1, and 0.016 x Mbh for set 2, respectively. Values 
at r sg (S and H) are used as the inner boundary conditions for 
the self-gravitating part. 

The emergent spectrum from r < r sg (for set 2) is shown in 
Fig-E]by a solid line (right), while relevant physical quantities 
are shown in Figure|5] For a presentation purpose, a simple T e g 
profile (i.e. T e g oc r~ i/4 ) is also drawn as a dotted line. 
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Fig. 4. Mid-IR to X-ray spectral modeling of Ton S 180. Open squares and thick striped lines are the observed data and the 
power-law fits described in Turner et al. (2002). Two solid lines are the spectral components of the dusty torus (left) and the 
non self-gravitating disc with a corona (right; Mbh = 10 68 M o , M = 1000LEdd/c 2 , and a = 0.002). Two dotted lines mean 
the spectrum (from r sg < r < r ou t) of the self-gravitating disc for y = -0.6 (where S oc r y ) and ff out = 0.02 (upper; one of 
the successful fits shown in Tabled) and that of a Keplerian, standard disc (lower). Total spectra of those three components 
are indicated by dashed lines: thick line is the spectrum with the self-gravitating disc mentioned above, while the case with a 
standard disc is indicated by a lower curve. The upper dashed line exhibits a total spectrum with a different parameter set for the 
outer self-gravitating disc (y — 0, a out — 0.02). 



As to the broken power-law description for the X-ray 
data, different combinations of the photon index F and 
breaking energy are deduced for different observations. 
In order to assess how the different combinations affect 
the generation of the best sets of parameters, we have 
also performed the same least-square procedure with the 
X-ray data replaced by the broken power-law description 
in Vaughan et al. (2002). The following four parameter 
sets provide best fits; [log(M B H/M ), M/(L E dd/c 2 ), a] = 
(6.8, 1000, 0.002), (6.8, 1000, 0.004), (7.0, 562, 0.002), and 
(7.2, 178,0.002). The goodness of the fit to the data shows a 
slightly shallower distribution than the fit to the Turner et al. 
(2002) data. In the rest of this study, we emphasize the fit to 
the data in Turner et al. (2002), where the broadband data were 
obtained quasi simultaneously. 

3.2. Spectrum of the non-Keplerian, self-gravitating 
disc 

In the standard theory of radiatively cooled accretion discs 
(e.g., Pringle 1981), the effective temperature T e ff (close to the 
surface temperature) is given by 



0-7WV) = v t S r— , (11) 



where Q is the rotation frequency. Again, the turbulent vis- 
cosity is chosen to be the standard prescription: v t = ac s H 
(Shakura & Sunyaev 1973). The sound speed c s is found from 
the requirement of hydrostatic equilibrium, namely 

c 2 = D^H 2 + AnCLH. (12) 

At r > r sg , the disc becomes self-gravitating, namely Q 
is expected to depart from its Keplerian value Qk- In order to 
compute O, we assume that the rotation curve has only a grav- 
itational origin, that is 

rrf(r) = ^S-rf-(r), (13) 

where g* sc is the contribution due to the whole disc, including 
the inner slim disc. In general, at large distances, we expect 
g^ lsc < with the consequence that Q. > Qk- This means a 
possible reduction of the gradient, \dQ./dr\, and subsequently a 
change of the effective temperature [eq. il 1H with respect to 
a standard disc. It is worth noting that a similar idea has been 
proposed by Lodato & Bertin (2001) in order to explain the 
IR-excess of T-Tauri stars. 

Since the self-gravitating disc should be gravitationally un- 
stable, other sources of heating, e.g., heating by gravitational 
instabilities (Adams, Ruden & Shu 1988; Lodato & Bertin 
2001), may overcome the heating by turbulent viscosity (see 
Lodato & Rice 2003). To take into account this possibility, we 
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Table 2. Successful solutions for the self-gravitating disc. 



r 


Q'out 


M^O < r ou t)/ Mbh 


0.3 


0.002 a 


7.0 


-0.6 


0.02 


1.4 


-1.5 


0.1 


0.4 



Namely, constant a in the entire disc. 



allow the parameter a in the self-gravitating part to increase 
linearly with r, from 0.002 (the value used for the inner disc in 
§ 3. 1) at r sg to a out at the outermost radius of the self-gravitating 
disc r out . 

Since there is no simple formula for radial self-gravity, we 
have determined g* sc numerically using the accurate Poisson 
3D-solver, which is described in Hure (2003) for discs with 
various shapes, sizes and surface density profiles. This has been 
performed assuming that H oc and E oc R? , where /3 and y are 
input parameters. Here, H and E are continuous at the transition 
between the two regimes of the disc. The radial extension of 
the self-gravitating disc r out and the viscosity parameter a at 
that radius a oat are also input parameters. Those values (J3, y, 
r out , and o- ut) are determined so that the total spectra (§3.4) 
provide acceptable fits to the observed one. 

For each model, the effective temperature T^{r) is derived 
from the resultant rotation frequency using Eq. dl It . The spec- 
trum of the self-gravitating disc is then computed in the same 
way as the inner slim disc (i.e., with the effect of electron 
scattering, even though the effect is small in the outer self- 
gravitating part). 

When we add the spectra of the outer, self-gravitating disc 
onto the inner disc spectrum, it is found that the combined spec- 
tra based on set 2 (Table0 generally provide a better fit to the 
observed optical/UV spectrum than those based on set 1 . Thus, 
we hereafter discuss the models obtained for set 2 alone as an 
example. 

For set 2, three possible combinations of parameters are 
shown in Table |2] The best parameters for ft and r out are com- 
mon to them: (3 - 1 and r out = 3 x 10 4 rs c h (= 10 x r sg ). For 
one case (y = -0.6 and a oal = 0.02), the spectra from the self- 
gravitating part (upper dotted line) and the total spectrum (thick 
dashed line) are presented in Fig.@] Those three combinations 
exhibit almost identical spectra (deviation of vL v is 0.05 dex 
at most). Other combinations are also acceptable. The cumu- 
lative disc mass Ma sc {r < '"out) is not negligible, ranging from 
0.4Mbh to 7.0Mbh- For the intermediate case, for instance, the 
total mass, M B h + M disc , is about 2.4M BH = 1O 7 2 M , which 
is consistent with the central mass estimated by Hfi and [O III] 
line widths. It indicates that Mbh -estimations using line widths 
overestimate the BH mass systematically. To assess how the 
spectrum changes with different y, the total spectrum obtained 
with y = 0, which is less preferable, is indicated by the upper 
dashed line. 

Figure[5]presents several physical quantities of the disc for 
the three successful parameter sets (Table |2} as a function of 
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Fig. 5. Variation with radius of some disc parameters (the in- 
ner, non self-gravitating part is based on set 2); half surface 
density E, rotational frequency normalized to the Keplerian 
value CI/Qk, effective temperature T e g, and Toomre Q- 
parameter. The solid lines are computed with y = -0.6 and 
ffout = 0.02, while dot-dashed and dashed ones are computed 
for y = 0.3 with a constant a and y = -1.5 with ff out =0.1, 
respectively. Two accretion regimes are detached at r = r sg (= 
3000r Sch ). 

radius; E, O (Ea. ll3t and r e ff. Vertical optical thickness of the 
flow, t [t = E x (/c a bs + K es ) > E K es = 0.4 (E/g cm -2 )] is much 
larger than unity everywhere. For the case with a very massive 
disc (y - 0.3), Q. becomes flat and even starts to increase with 
radius due to the strong radial self-gravity. This produces a lo- 
cal decrease in T e ff [oc (dCl/dr) 05 ; eq. dl 1H at a small range of 
radius. However, the temperature profiles T e ff(r) for the three 
cases are mostly the same, and thus they produce almost iden- 
tical spectra. 

3.3. Emission from the dusty torus 

The dusty torus component is supposed to be responsible for 
the Mid-IR (~ 30yimi) to ~ 2{im emission. In general, the Near- 
IR flux from AGN varies with time following the optical flux 
with a time lag of about a month-year (e.g., Clavel, Wamsteker 
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& Glass 1989; Glass 1992; Nelson 1996). This implies that the 
dusty torus is located quite far from the BH, and that it is heated 
by the radiation from the vicinity of the central BH. 

A detailed modeling of the emission of the torus is beyond 
the scope of this paper. Here, we adopt a simple power-law 
description with a cut-off. Theoretical arguments on the dust 
sublimation as well as observational data (as a spectral bump 
around 3 fim) indicate that the maximum temperature of the 
torus is ~ 1000- 1500K (e.g. Efstathiou, Hough & Young 1995; 
Kobayashi et al. 1993; Pier & Krolik 1993). It corresponds to 
wavelength of ~ 2.4/rni [i.e. hc/(4kT) with T = 1500K], or 
v 10 141 Hz (i.e. 4kT/h). 

We assume that the fluxes at 25 fim and 12/im come from 
the dusty torus alone, while the near-IR flux is partly due to the 
self-gravitating disc. Then, the simple model we apply is 



VLy = 10 



38.45, ,0.45 



exp(-v/10 141 Hz)erg/s. 



(14) 



We are confident that this is neither the only nor the best possi- 
ble description. The spectral index inferred for this object dif- 
fers from the typical mid-IR-near-IR index of bright AGNs 
(vL v oc y~ - 4 ; Neugebauer et al. 1987; Polletta & Couvoisier 
1999; see also Grupe et al. 1998). 

We can specify the location r t of the inner torus (see FigQ 
by considering that the torus is in thermal balance (heating by 
radiation from the inner disc is balanced by radiative cooling): 



r, 5.3 x 10 1 



L 



0.5 



10 45 erg/s/ \1500K 



-2 



cm. 



(15) 



where L and T are the bolometric luminosity of the central 
disc, and the gas temperature at the innermost radius of the 
torus, respectively. For the set 2, we have r t 5.3 x 10 17 cm(= 
0.17pc) =2 2.8 x 10 5 rs c h, which is about 10 times larger than the 
outer radius of the self-gravitating disc r out . Thus, gravitational 
attraction on the self-gravitating disc by the torus should not be 
crucial, unless the torus mass exceeds the BH mass by some 
order of magnitude. 

3.4. Mid-IR - X-ray spectral energy distribution 

The thick dashed line in Figure |4] shows the broad-band, mid- 
IR to X-ray spectral modeling for Ton S 180, comprised of the 
three components described in the subsections above. We see 
that the optical component can be interpreted as the thermal 
emission from the self-gravitating, non-Keplerian disc. The 
models for the self-gravitating disc with y and ar out listed in 
Table [5] provide equally good fits to the observed spectrum. 
The upper dashed line (for y — and a oM = 0.02) is drawn 
to show to what extent a different value of y changes the total 
spectrum. 

If the outer self-gravitating disc is replaced by a Keplerian, 
standard disc, the resultant total spectrum (lower dashed curve) 
is not adequate to explain the observed optical/near-IR spec- 
trum. In other words, flux enhancement at the outer region is 
necessary. This enhancement depends on the resultant rotation 
law, as well as a(r), at the outer self-gravitating disc. 

We demonstrate that a broadband spectrum of one NLS1, 
Ton S 180, is indeed well fitted by a summation of an inner, 



non self-gravitating disc, an outer, self-gravitating disc and a 
dusty torus. Detailed mid-near-IR observations, including tem- 
poral studies (e.g. Glass 1992; Nelson 1996), will provide a 
more accurate model of the dusty torus. Subtraction of the torus 
component from the total (disc and torus) spectrum will en- 
able us to derive stronger constraints on the size of the self- 
gravitating part, and mass of the whole disc, etc. High spatial 
resolution obtained with ground-based, mid- to near-IR inter- 
ferometers with the large telescopes, such as OHANA (Mariotti 
et al. 1996; Perrin et al. 2000) and VLTI, will contribute to solv- 
ing these problems. 

4. Implications and limitations of the present 
study 

As is obtained in §3.2, the integrated disc mass is about 
(0.4 - 7.0)M BH (Tab. Although the Toomre Q parameter 
[= Q.c s /(nI.G)] remains greater than one in the whole disc 
(Fig |5j, it is quite close to unity at outer radius for larger y. 
This fact supports the introduction of the additional heating 
due to gravitational instabilities (e.g., increasing a used here; 
see Lodato & Rice 2003). Consensus on the amount of such 
heating has not yet been reached. Further understanding on this 
topic, which is related to a out in this study, will allow us to con- 
strain y (where S oc r y ) and the disc mass more severely. 

The size of the broad-line region (BLR), /?blr, for 
Ton S 180 is expected to be of the order of 100 It. d., us- 
ing the empirical relation between the optical luminosity and 
^blr (Kaspi et al. 2000). This radius corresponds to 0.085 pc 
= 2.6xlO I7 cm;equivalently, 1.4xl0 5 r Sch forM BH = 10 6 8 M o . 
Thus, the broad-line region is likely located between the in- 
ner radius of the torus (r, ^ 3 x 10 5 rs c h) and the outer radius 
of the self-gravitating disc (r out = 3 x 10 4 rs c h)- However, the 
constraint on r out derived in § 3.2 is rather weak. If the outer 
massive disc extends at r » r out , the clouds in the broad-line 
region may not be in Keplerian motion, which will lead to an 
uncertainty for the Mbh -estimation with H/3 width (e.g., Krolik 
2001). 

Since the main purpose of this paper is to report the discov- 
ery of the spectral window for observing a self-gravitating disc, 
and to demonstrate a spectral fit to the observed spectrum of a 
super-Eddington accreting AGN, we did not treat the following 
two issues very accurately: (i) the vertical self-gravity (the last 
term in the right-hand-side of Eq. I12> in the inner, non self- 
gravitating disc, and (ii) the effect of realistic absorption opac- 
ity in deriving Pmid- P iane (§2.1). Also, possible outflow from the 
disc is not taken into account either in the inner slim part or in 
the outer self-gravitating part. Outflow/evaporation at the outer 
region may contribute to the formation of BLR clouds or broad 
absorption line (BAL) clouds. These issues will be examined 
in the future. 

5. Conclusions 

We have presented the outermost radius of the non self- 
gravitating accretion disc (r sg ) around a super-massive black 
hole (BH), i.e. the inner edge of the self-gravitating disc, as a 
function of the BH mass Mbh, accretion rate M and viscosity 
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parameter a. We then showed the corresponding wavelength 
(A sg ) that is relevant to the emission from r sg . 

When M < 70 L EAd /c 2 (for M BH = 1O 7 M ), A, g is roughly 
proportional to M~ as for fixed M BH and a. With a higher M, 
A sg reaches and stays at ~ 4000 A for a = 0.1. Interestingly, A sg 
is as small as ~ 1500A for the lowest a case in this study (i.e. 
a = 0.001). Therefore, the continuum emission from ~ 2/j.m to 
/l sg arises from the self-gravitating part of the super-Eddington 
accreting disc. Thus, we have discovered a spectral window for 
observing the self-gravitating disc, which has only been studied 
by maser spots at the parsec-scale for a few nearby objects. 

Next, we demonstrated that the mid-IR to X-ray spec- 
trum of a bright, well-studied Narrow-Line Seyfert 1 galaxy, 
Ton S 1 80, is indeed well fitted by the spectrum arising from the 
following three components: an inner slim disc (with a corona), 
an outer, self-gravitating non-Keplerian disc and a dusty torus. 
Comparing the model spectra of the slim disc with the ob- 
served UV-X-ray one, the following parameters are favored: 
M BH ~ 1O 68 M , M * 1000L E dd/c 2 , and a « 0.002. In the 
model for the outer, self-gravitating disc, we allow the viscos- 
ity parameter a to increase with radius, and an outer radius of 
3 x 10 4 rsch is inferred. Various profiles for the surface density 
X are acceptable, e.g., from X oc r° 3 to £ oc r~ L5 , depending on 
a at the outermost radius of the disc. 

The accretion disc is quite massive, and the disc mass can 
be (0.4-7.0)Mbh- The total mass, BH mass plus the entire disc 
mass, is found to be about 2.4 M BH = 1O 72 M if E oc r~ oe \ 
which is consistent with the central mass estimated by H/3 and 
[O III] line widths. This indicates that MfiH-estimations using 
line widths systematically overestimate the BH mass. Although 
the disc mass is quite high relative to the BH mass, the Toomre 
Q parameter remains greater than one in the whole disc, so that 
the disc is marginally stable. 

Moreover, r sg can be less than the photon trapping radius 
for M > 10 33 LEdd/c 2 . Thus, a self-gravitating, optic ally- thick, 
advection-dominated accretion disc is expected to appear in the 
extremely high accretion rate regime. 
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